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RECORDING SYSTEM CLASSIFICATION 


Analog Recording Systems — CenNtTinuvovS 


- Information (message) signal to be recorded has infinite number of 
amplitude levels that change continuously with time. | 


- Typical Requirements: high signal-to-noise ratio (SNR), low dis- 
tortion, and low cost. 


Digital Recording Systems — Disceem LEMS 


- Information (message) signal to be recorded has finite number 
(usually 2) of amplitude levels that change at discrete time points. 


- Typical Requirements: high reliability (low probability of error), Tast 
access {to recorded information, and low cost (S/Mbyte). 
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RECORDING SYSTEM APPLICATIONS 


Enclosure 


Disk Stack 


Electromagnetic 
Actuator 


- Rigid Disk Drives 


Swing Arms 


- Rigid disk file components. ( Rem Mee & Daniele) 


Jacket 


- Flexible Disk Drives 


Stepping Motor ; 
Flexible diskette and head accessing system. (From Mee & Dan | he) 
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RECORDING SYSTEM APPLICATIONS (contd.) 


Audio Recording Systems 


- Analog Audio Recording 


Cassette 


Record/Replay Heads 
Capstan and 


Pressure Roller Sx) 


Audio cassette recorder. (Rom Met ePem 


- Digital Audio Recording 


RECORDING SYSTEM APPLICATIONS (contd.) 
Image Recording Systems 


- FM Video 


Rotating-Head Drum 


Read-Write Head 


Cassette 


on Helical-scan video ; te) 
recorder. ( Fanrw Me2e © Dames 


- Digital Video 


RECORDING SYSTEM APPLICATIONS (contd.) 


Instrumentation Systems 


Magneto-optical Recording Systems 


Signal 
Photocetectors 


Focus and TracKing 
Detectors 


Tracking Mirror 


Megnetoontical 


reed-write need for recora- 


Components for an optical | 
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DIGITAL MAGNETIC RECORDING SYSTEMS 


General Requirements 
- Low Probability of Error 


- Fast Access to Recorded Data Reve Time SIGNAL PLocEssiA 


- Low Cost 


Parameters of Interest 


- Linear Density, number of bits per unit length along a track 
(bits/inch) 


- Track Density, number of tracks per unit length (tracks/inch) 


- Areal Density, number of bits per unit surface area (bits/inch’); 
product of track and linear density. 


- Volumetric Density, number of bytes per unit volume (MBytes/ cu. 
ft.) 


STATE-OF-THE-ART IN TPI AND BPI 
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Track Density (t/mm) 

(a) Areal Gensity progress for magnetic recording 

oa (bd) re tionships of linear and track densities. Tne 
nslty curves for eudio and video recording are obtained from 
the approximation thet the shortest recorded wavelength 18 
uivalent to two recorded bits. HV, home video tape; PY, 
ofessional! video tape; RD. rigid disk; DT, date tepe; FD, flex- 
Die sie AT, audio tape (DAT, digital audio tape); MO, 
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Gaussian distribution 


If X represents the sum of a large number of independent random components, and if each | . 
component makes only a small contribution to the sum, then 
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Rayleigh distribution 
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here X and Y are independent gaussian r.v.s with zero mean and 
variance c’, then 
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TABLE T.6 
GAUSSIAN PROBABILITIES 


The probability that a gaussian random variable with mean m and variance a? will have 
an observed value greater than m + ko is given by the function 
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called the area under the gaussian tail. Thus 
P(X > m + ko) = P(X < m— ko) = Q(k) 
P(|X — m| > ko) = 2Q(k) 
Pim <X <m+ko) = Pm—ko <X <m)=¥,- Oth) 
P(|X — m| < ko) = 1 — 2Q(k) 
P(m—kyo < X <mt+k, a) =1— Q(k,) — Q(z) 


Other functions related to Q(k) are as follows: 
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All of the foregoing relations are fork > O. If k < 0, then 
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Numerical values of Q(k) are plotted below for 0 < k < 7.0. For larger values of k, Q(k) 


may be approximated by 
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which is quite accurate for k > 3. 
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Operation 


RESPONSE OF LINEAR SYSTEM 


PROBLEM OF DETECTION IN MAGNETIC RECORDING 


PULSE RESPONSE AS A FUNCTION OF LINEAR DENSITY 
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INTERSYMBOL INTERFERENCE (ISI) IS A KEY EFFECT OF 


HIGH-DENSITY RECORDING 


increased isi <== 3 lowered signal-to-noise ratio (SNR) 


lower SNR s===—> higher error rate 


SIGNAL PROCESSING AND CODING METHODS 


- Peak Detection Method 

- Partial Response Methods 

- Equalization Methods 

- Maximum-likelihood detection (Optimum) 
- Codes for Peak Detection 


- Codes for Partial Response Methods 


- Signal space coding (trellis coding) 
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FOURIER TRANSFORMS 


Definitions 
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Transforms 
Function Xt) V(f) 
Rectangular I] (:) tT sinc ft 
Triangular a(*) t sinc? ft 
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for x from 0 


sinc? x 
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TABLE T.5 
PROBABILITY FUNCTIONS 


Binomial distribution 


Let the discrete r.v. J be the number of times an event A occurs in n independent trials. If 
P(A) = a, then | 


j = na c? = na(l — a) 
Ifn > 1,a<«< 1,and m = na remains finite, then 


Pi) = e7™m'/i! 


Poisson distribution 


Let the discrete r.v. J be the number of times an event A occurs in time T. If P(A) = 
u AT «1 ina small interval AT, and if multiple occurrences are statistically independent, 
then 


Pi) =e (uTVi! i= pT 


Uniform distribution 


If the continuous r.v. X is equally likely to be observed anywhere in a finite range, and 
nowhere else, then 


Sinusoidal distribution 
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Least-squares 
storage-channel — 
identification 


Pulse (dibit) and step (transition) responses for 
magnetic-storage channels are important for 
detection-circuitry design and for comparison of 
various media, heads, and other channel 
components. This paper presents a least- 
squares procedure that can be used to identify 
the dibit and transition responses from 
measurements of the read-head response to any 
known data sequence written on the medium. 
The method yields significantly higher-quality 
estimates for the dibit and step shapes than 
does determining these same characteristics by 
measuring the average response to isolated 
transition or by performing a Discrete Fourier 
Transform (DFT) on the response to a 
pseudorandom data pattern. The new method 
can be implemented off line but also can be 
made sufficiently efficient to be implemented 
with a microprocessor for use in self-optimizing 
(adaptive) channel detection circuitry. 


1. Introduction 

Storage-channel identification 1s the measurement and/or 
computation of the characteristics of the read-back channel 
in a data storage device, such as a magnetic disk, magnetic 
tape, or optical disk. The identified characteristics are most 
often the channel’s response to a step input (the “transition” 
response) or to a pulse (the “dibit” response). These 
characteristics are important for many purposes, such as the 
design of the detection circuitry (especially for equalizers and 
©Copyright 1986 by International Business Machines Corporation. 
Copving in printed form for private use is permitted without 
payment of royalty provided that (1) each reproduction is done 
without alteration and (2) the Journa/ reference and IBM copynght 
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for maximum-likelihood detectors), for determining the 
maximum data density of the device, and for comparing 
various media, heads, and other channel components. 

This paper presents a least-squares procedure for 
identification of the linear time-invanant filter that most 
closely approximates the desired step or pulse responses. The 
storage device is excited with a known data sequence, and, 
later, the read-head response to the known sequence is 
measured (or digitized) at regular intervals. The resulting 
measurements are then processed via the least-squares 
procedure to determine the step and/or pulse responses. 

The resultant estimates of these responses are of 
significantly higher resolution (higher quality) than those 
produced by previous procedures, such as measuring the 
average response to isolated transitions (or isolated dibits) or 
computing the Discrete Fournier Transform (DFT) of the 
response to some known (usually pseudorandom) data 
pattern. Furthermore, the new method, although based on a 
linear model of the channel as presented here, can indicate 
the average accuracy of the linear model over any data 
pattern, thus indicating the presence of potential 
nonlinearities in the responses, unlike the aforementioned 
methods. The degree of agreement between the linear model 
and measurements can be useful in determining the data 
rates at which various data detection methods do and do not 
apply. 

Section 2 defines in more detail the quantities used in 
channel identification and the least-squares procedure, and it 
compares the quality of estimates of the new and previous 
procedures. Section 3 studies some details of the solution 
and displays the results of the new procedure for several 
measurements taken from actual storage devices, including 
magnetic disks with thin-film heads, tape systems with 
magnetoresistive heads. and optical disks. Section 4 is a bnef 
conclusion. Appendix A extends the channel identification 
procedure to apply at any digitizer sampling rate (an integer 
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atio of the sampling to data rates is assumed in the main 
body of the paper). Appendix B discusses streamlining of the 
:2ast-squares procedure for possible use with adaptive 


Jctection methods, while Appendix C discusses the detection 
of nonlineanities. 


a Storage-channel identification methods 
“Phis section mathematically defines and analyzes the 
4uantities and procedures used in storage-channel 
tdentification. Figures 1(a) and 1(b) summarize the 
definitions used throughout this section. 


e lariable definitions 
1e@ read-back channel and associated identification 
-Tameters are illustrated in Figures 1(a) and 1(b). The 


:ntinuous read-head output signal, a(t), can be modeled in 
ae of two ways [1]: 


dig = ¥ x,A(t — kT) + u(t), (1a) 


AC. H= ¥ s,4(t -— kT) + u(t), (1b) 


were A(t) and h(t) are the unknown linear time-invariant 
pulse and step responses, respectively, and u(t) denotes an 
wacorrelated, additive, zero-mean noise,t x, takes on the 


~~ tajues +1 (or +1 and 0 for some optical storage systems), 


otresponding to 1’s and 0’s, respectively, in the stored data 
sequence at time KT, 1/T is the data rate, and k is an integer. 
!3 Equation (1b), s, can take on the values +2 or 0 (+1 or 0 
fer optical) according to the relation 


SF XT Xa (2) 
Likewise, one determines for a linear channel 
hit) = A(t) — A(t — T). (3) 


[tis a property of the method presented that the estimates 
also obey Equation (3); however, it is sometimes informative 
to separately identify A(1) and h(t), rather than identify only 
one and compute the other from it. It 1s assumed that d(r) is 
digitized at some rate 7.,, such that 


en (4) 


where p is an integer (=1) oversampling factor. This 
restriction is relaxed to a rational fraction in Appendix A. 
The sampled read-head output is then, with f= mT, in (1), 


d(mT,) = ¥ x,A(mT, — kT) + u(mT,) 
k 


lI 


>» x,Aftm — kp)T,) + u(mT,) (Sa) 


k 


ven though the assumption that the noise 1s additive may not be completely true in 
-acuice. our objective 1s to find the values for the parameters in such a model that 
most closely approximate the measured responses, and deviations from such a model 
appear in the final results of the method in this paper. 
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Summary of storage quantity definitions (a) for pulse responses and 
(b) for step responses. 


or? 


dmT,) = ¥ s,h(mT, — kT) + u(mT,) 


k 


=) s,h[(m — kp)T,] + u(mT,). (5b) 


The channel is estimated by 


amT,) & ¥ x,w(mT, — kT), (6) 
k 


where w(f) is a linear filter response whose sampled values at 
times mT, are to be computed via the channel identification 
procedure [ideally w(t) = A(1)]. Likewise, for the step 
response, the estimate is 


d(mT,) & ¥ s,wmT, — kpT,). (7) 
kK 


We also define an error signal 
emT,) & d(mT,) - d( mT,,). (8) 


As an example, note that, if x, or s, 1S a sequence 
corresponding to an isolated pulse or transition input, (Sa) 
and (5b) reduce to 


aAmT.,) = h(mT,) + u(mT,) (9a) 
or 
d(mT,) = h(mT,) + u(mT,), (9b) 


respectively, the desired pulse shapes in noise. Then, (777) 


t The reader may note that (Sa) and (5b) are equivalent to p subchannels, each at 
spacing J. this observation is exploited to reduce computation in the new procedure 
in Section 3. 
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and w(mT_,) can be estimated by the averages 


wOmT)) = - XY dmT,,k), (10a) 


A=] 
] n 

wm) = ; Y d(mT,k). (10b) 
k= 


where the index k denotes the kth experiment. That is, one 
measures the response n times and averages, which is the 
basis for the aforementioned isolated step and dibit 
identification methods. Some deficiencies of the estimates 
identified via such isolated step or pulse methods are 
discussed later. Equations (9) and (10) were given only to 
verify the utility of the definitions in (1)-(8). We now 
proceed with a discussion of the least-squares channel- 
identification procedure. 


e The application of least squares 

In the least-squares identification procedure, a known data 
pattern is written on the storage device. The w(mT_,) are 
chosen to minimize 


l 
t= D my7,), (11) 


m=] 


where e(mT,) is given in (8). If we denote W,,, by the M x | 
column vector 


4'(0) 
Weal  s _ |, (12) 
wi{(M — 1)7,] 


then the solution to (11) 1s conveniently written [2] 


_ in cal 
Wa ( 2 Kink) ( > Xu md mT,)), (13) 


m=] m= | 


where ’ denotes transpose, and 


Xm 
X atm 2 ° (14) 
Xm-M+) | 


for p = 1. There are p — | zeros between entries 1n (14) if 
p> |. We have further assumed that A 1s large enough to 
span the nonzero extent of the pulse (step) response in 
intervals of sampling periods or M7, = NT data penods 
containing p samples each, M = Np. Equation (13) can be 
rewritten 


Wa = Rie M1? (15) 
where 
| / 
Ry = | y) Xm *X Mom ’ 
m=] 
| i 
Prag a) » Xap mA MT 4). (16) 
m=} 


A similar expression holds for the step response. with x's 
replaced by s’s and \"s replaced by w.,’s in the solution. Note 


teen 
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that M x M matnix inversion is explicit in (13); however, 
because of the special structure in this problem, no matrix 
need ever be inverted directly. For more details, see Section 
3 and especially [2]. 


e A performance measure 
The mean of W,,,can be easily determined as 
h(0) 

EW) = Ay = , — (17) 

h{((M — 1)T,] 
the desired solution, when the above least-squares method is 
used. The Norm Tap Deviation is a mean-square measure of 
statistically how far the estimated W,,,1s from H,,, and is 


also easily computed, if u(t) 1s white (spectrally flat over the 
frequency range of interest), as 


I 


Bay) = Ell Way) - Hy] ae trace (Ryo, , (18) 
where 
o, & E[u(KT,)’]. (19) 


We show in the next few sections that both the isolated 
transition (or dibit) and DFT methods are special cases of 
the general least-squares method with very special 
restrictions on the input sequence and on M and /. Thus, we 
are able to use (18) as a performance indicator for those 
methods as well. 


e /solated transition example and analysis of resolution 
AS an example, once again consider an isolated dibit; then 
X sm has only one nonnegative entry per column and (13) 
reduces to (using generalized inverses, see [3]) 


* M 7 
Way = . : (20a) 
a T,) . 


A string of n “isolated” (far enough apart) dibits occurring 
within a large data record (length /) has a least-squares 
solution, 


d(kMT, + T,) 


1 n-| 


(20b) 


Ws = hp ° | 

k=0 | AKMT, + MT,) 
that is exactly the same as the isolated pulse solution in 
(10a). The least-squares identification procedure is more 
general in that the input need not be an isolated transition or 
dibit. 

Equation (18) allows us to compare the quality of the 
least-squares estimates of H,,,(W,,,) for different input 
sequences. Note that, for a string of 7m isolated (MT, apart, 
so / = Afn) inputs. one determines for white zero-mean wu 


. ’ > M 2 
EW Wy) — Ay = Son (21) 
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Pseudorandom sequences are generally desirable [4, 5] for 
channel inputs because of their broadband spectral response. 
An identity for R,,,can easily be determined, if / = length of 
the pseudorandom sequence,§ as (see [6-8]) 


Ry = + Wag — Lelia, (22a) 


where 1,, 1s an Mf x | vector of M ones. One can also easily 
show that 


he av tt 1 
Rut = Ty (nt 5s , hala) (22b) 


Thus, (18) becomes, for a pseudorandom sequence of length 


‘f repeated n times, 


- oh (ee 2,1 92 03) 
“aM +l (n- DM 4 fn 


-orn = 1, there is an improvement of (/ + 1)/2 with 
“espect to (21). As 1 increases to a large value, there is an 
‘mprovement by a factor of M in estimate quality, or 
‘quivalently, 44 more digitized outputs from isolated dibits 
“ust be processed in the isolated dibit identification schemes 
:@ get the same resolution estimates as those produced by 
2ast squares with a pseudorandom length-/M input. For 
oversampling (p> 1), the companson favors the 
“seudorandom input by the same amount. Heunistically, 
“hen using pseudorandom or “scrambled” data in channel 
stentification, the input 1s more spectrally “rch” and all 
‘yequencies are more equally weighted than when a single 
sulse is used. The resulting flat nature of the spectrum 

“osults in the inverse autocorrelation matnx being close to an 
“Jentity which makes 6,,,in (18) smaller (better). When x, 
has a flat spectrum, s, does not have a flat spectrum, but a 
similar slightly more complex argument can be given to 


justify the least-squares improvements. 


In practice, it may not be difficult to average the extra data 
for the isolated input method. However, there is another 
very practical advantage of using more random data, as was 
first noted by C. M. Melas [9]. This is that in the isolated 
transition or isolated dibit methods, the AGC (Automatic 
Gain Control) must be removed from the channel to prevent 
the sudden change in energy associated with the isolated 
input from suddenly varying the gain parameter of the AGC. 
Then, the identified pulse charactenstics will not include the 
effect of the AGC. This effect can commonly be more than a 
simple gain factor and is determined by the bandwidth and 
tracking rate of the AGC. 


© Comparison with frequency-domain methods 
Another more recent method used in storage-channel 
identification is [4, 5, 10] to compute the DFT of the 


bes ee wa 
> f.en when the output is oversampled. we show later that the only autocorrelation 
Matria of interest 1s at the data rate; thus all of the analysis here 1s also valid for p> 1. 
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response to some prescribed pattern written on the media. In 
order to invert the DFT to get a time-domain estimate of the 
pulse response, one must first divide the measured DFT by 
the DFT, including phase, of the input before the inverse 
DFT, which [10] also observes. Using this last restriction, 
one can also generalize the methods of (4, 5, 10] to estimate 
the channel response for any inputs, including the +2, 0 
normally associated with identification of the step 
(transition) response. 

Nevertheless, with the division by input spectra. the 
frequency-domain method is the same as the time-domain 
least-squares method of this paper if 44 = /, and as we shall 
see, the case M = / gives very poor estimate quality. In the | 
case that u(mT_,) is white and Gaussian, the least-squares 
method (see [11]) achieves the famed Cramer-Rao bound for 
a fixed / and M; that is, no other estimator has higher 
resolution for the given data. If the assumption on u(f) is 
just white (not also necessarily Gaussian), then the least- 
squares estimator 1s a Best Linear Unbiased Estimator 
(BLUE) [3]. 

Theoretically, the difference between the DFT technique 
and the time-domain least-squares method can be quantified 
via the following analysis. It is usually wise to pick / </so 
as to introduce more noise averaging, or equivalently, to 
make the Cramer-Rao bound lower for fewer parameters. 
Generally speaking, in any estimation scheme. we desire 
1 > M to get good quality estimates. Nevertheless, picking 4 
too small can introduce extraneous harmonic distortion in 
the estimated step response. The time-domain least-squares 
method can be rewritten as that W’,, , that minimizes [2] 


En = Extasy = | aul | (24) 
where 
Eng = Fy — Xap a4 (25) 
and 

A kT ,) x, 
a), = ; Xi = 7 ds (26a) 


d{(k -1+1)T,] 


Na+ 


where p — | zeros can be inserted between nonzero entries in 
X,, and 


X pin = (Xie Xie 07% Xitel: (26b) 


The DFT-based method is a special case of a linear Af x / 
transformation on e,,,. that is, let 


Evra = Plats» (27) 


mene 4 


where ¢ is an Mf X / (possibly complex) matrix representing 
the linear transformation. Then 


* — ~_* * 5 
E Waa = 7 PPE v7) 5 (28) 


where * denotes conjugate transpose. If ¢ is a unitary 
transformation (¢*@ = /). then 
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periods. 


eo 


NEw il = lea’, | (29) 


and the minimized ¢,,,1s obtained by 
Ey) = CE, - (30) 


In the DFT methods of [4, 5], the matrix ¢ is chosen, under 
the very special assumptions that M = land the input is 
periodic ( pseudorandom) of length |= M, as 


l l 


1 1 gor Pies ual 


VM 1 Jw, T J wag (M—-1)T 
e - @ 


Oa i=0,---,M- 1. (32). 


¢,, can easily be shown to be unitary [12]. so the relation in 
(29) holds, apparently yielding the time-domain least-squares 
solution. #%, 1s the inverse DFT in this case. However, in the 
time-domain method of this paper. .\f1s much less than 

/ to average the effects of noise and other nonideal effects. 
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Using our performance measure in (18) and (23) (n = 1, 
| = M) again, one determines the estimate quality as 


o., = +1 o, ’ | (33) 


while the general formula for a pseudorandom sequence of 
length / with M parameters is 


2M+MI-M , 
G+ il=M+i 


O47) = 


(34) 


Substitution of / = 10M, a good practical rule of thumb, into 
(34) yields the advantage 


(0.9/ + 1) 


A. 210.91 +1) _ 
(0.9/ + 2) - 


ae re ee 2 


6,,, 0.097 + 0.2/— (35) 


Even for / = 1000, another reasonable number, the 
improvement in (35) is close to its limiting value of 20. This 
large improvement is typically evident when comparing the 
spectra of a pulse produced by the time-domain least squares 
and by the DFT method, as we have illustrated in Figures 
2(a) and 2(b). Note from the level of “frequency npple” in 
the DFT plot that the time-domain least squares is at least 
an order of magnitude improvement. Also note the lower 
“noise level” at higher frequencies with the least-squares 
identification procedure. It is also important to note that 

/= M =2'— 1 (ia positive integer) for a pseudorandom 
input, which, at least, requires special attention for efficient 
DFT implementation [12-14]. The reason for the two 
different lengths (A/’s) in Figures 3(a) and 3(b) is discussed 
later. 


e An averaged DFT identification scheme 

Here, we propose an averaged DFT method for the special 
case that / = nM, where n is an integer greater than 1, and 
the input sequence is periodic with period M. [The case of 
oversampling (p > 1) is identical for each of the subchannels 
(see Section 3).] There 1s a very special set of circumstances 
when the inverted matrix in (13) is Toeplitz and DFTs can 
be used. Generally, (13) is not Toeplitz and DFTs are not 
appropnate. This method 1s equivalent to least squares, as 
can be seen from the following. Define ¢, by 


Psy 0 - 0 
Oy a are (36) 
0 0 .: Dry 
Multiplication by ¢ is equivalent to n M-point DFTs 
performed on the n groups of Af inputs. Note that ¢, is 
unitary, 
bo; = I. (37) 


The least-squares estimates in the frequency domain are 
given for each frequency bin by 
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d 4A. f)x*(k) ; d 4(k, 7) 
p(k) —_ ik ee — —___ 
Same. A 


i=] 


k=0,---,M-—1, (38) 


where »(A) and x(k) are the 4/-point DFTs of W,,, and 
Vir ay- Fespectively. 6(A, 7) is the M-point DFT of the time 
senes d, in the ith of the ” groups. Equation (38) is really the 
average of n uses of the original DFT method, when a 
nenod-.\/ input is recycled to fill / time periods. Then, some 
averaging will be introduced, in the optimal least-squares 
.ense, into the DFT identification scheme. The method of 
36) and (38). because of (37), is equivalent to an /-point 
‘oast-squares time-domain procedure. Of course, an inverse 
FT on the quantities in (38) must be performed to obtain 
‘ie desired time-domain parameters, W,, ,. This method 
-2quires the unnecessary imposition of an integer ratio 
-estricuon on / and m, which ts not required in the more 
zneral and straightforward time-domain least-squares 
lution (13). : 


- 4 note on maximum-likelihood detection schemes 

“he identified responses can be used in Maximum- 
_1kelihood Sequence Detection (MLSD) [15, 16]. In this 
case, the Mean Square Error (MSE) is a more useful estimate 
=f performance than (18). It is shown in [17] that (given a 
certain input sequence) 


USE = Efe(mT,)] = oy) (39) 
where y,,,1S given by 

; -1 
ya = Lm Xa Ray Xa (40) 


and ’ denotes transpose. One also can show (see [17]) that 
ONS Sail (41) 


thus. the worst (because the desired value is a”) MSE after 4 
measurements 1S 


MSE oon =U, (42) 


which is exactly the value given by a length-/ 
pseudorandom sequence. In fact, it 1s shown in [8] that 
choices for x, other than length-A/ pseudorandom sequences 
can yield MSE between 0 and o, after Mf data points, while 
sll maintaining good (low) E[|| Wy, — H,, MV). Thus the 
length-.\f pseudorandom sequence may not be the best 
training sequence if MLSD 1s used. Some data with statistics 
equivalent to what is expected in actual use would be the 
best choice for MLSD and other similar sampling detection 
schemes. 


© Signal-to-noise ratio estimation 
The SNR for the read-head response can be estimated (when 
Mi</) by 
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(= N) Wall? 
D- fury 


SNR = : (43) 


where || ’,, ,||°/p is the signal power for the binary input to 
the pulse response, and &,, ,// — N is the noise power. 
However, one must ensure that data measured at the read- 
head output have NOT BEEN AVERAGED before 
digitizing to ensure a meaningful estimate in (43). Also, as 
Howell [4] has noted, that distortion in the measuring 
devices, particularly the nonlinearities in the CRT sweep rate 
if a storage scope is used, can add appreciable noise not 
inherent in the actual storage channel. Of course, such 
contamination would leave (37) as a measure of the mean- 
square distortion in the measuring procedure, rather than 
the desired channel noise + media noise + modeling mean- 
square errors. Even if measurements are carefully taken, (43) 
is usually more indicative of the levels of nonlinear 
mismatch to the model and can therefore be very useful in a 
evaluating the potential success or failure of advanced _ 
detection schemes. | 


e Determination of M | 
We have previously assumed that the order M@ (number of cs 


identified parameters) was overestimated or known a priori. 


However, the best quality estimate for / data points is given 
by the so-called “Minimum Description Length” principle of 
[18], which jointly estimates M and the corresponding VW, , 
for /-points. The improvement in the general storage-channel 
identification problem is negligible if / = 10.47. It 1s 
interesting to understand just what happens if A/ is chosen 
too small. Suppose /A(AT,,) # 0 fork <0, k > M. Then 
u(kT.,) can be modeled as the sum of white noise and the 
distortion caused by the neglected terms in h. This second 
distortion term is just a linear filter acting on the 
pseudorandom pattern. When oversampled, the output of 
such a filter is the product of its transfer function and the 
transform of the oversampled pseudorandom pattern. The 
response of the oversampled pseudorandom pattern can 
easily be shown to be maximum at multiples of 1/7. thus 
explaining why choosing .M larger in Figure 2(b) than in 
Figure 2(a) caused the “harmonics” to disappear. Of course, 
picking M too large as in the DFT methods has a far more 
distorting effect on the output because of the lack of noise 
averaging. Generally speaking, conservative values for M/ and 
/ are 15 bit periods and / = 10M, respectively. 


e Summary 

In this section, we have introduced the least-squares channel- 
identification procedure, compared its performance with 

other commonly used procedures, and found the least- 

squares method superior in the quality of estimates that it 

produces. We now turn to implementation/programming of 

this new procedure. 315 
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3. Efficient implementation of the off-line least- 
squares identification procedure 

The time-domain least-squares solution 1s described using a 
matrix inverse in (13). This matnx can be large, requiring 
large storage and long processing time in an off-line 
computer program implementing the inversion. However, 
matrix inversion can be avoided to simplify the 
determination of W’,,,. This section describes several special 
features of the least-squares procedure that can be used to 
reduce considerably the computation and storage in an off- 
line implementation. Such simplifications could also become 
important if the characteristics of each particular storage 
device, and possibly at several different radii on each, were 
to be computed during the manufacturing process either for 
identifying defective devices or for optimization of the 
channel-detection circuitry for each particular unit. An 
efficient on-line procedure, similar to that of [8]. 1s suggested 


in Appendix B. 


e Subchannels 
In most cases of practical interest. the oversampling factor p 
is greater than one. Then, one wnites m7,= ni + 17, where 
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m 
r= (7) - 


where | denotes the “greatest integer less than.” and i takes 
the values 0, ---, p — 1. Equation (Sa) is rewntten [(5b) can 
be similarly rewritten] 


AnT + iT,) = ¥ x,h[(n — k)T + iT,]) + u(nT + iT,). — (45) 
k 


The index / has no effect upon the convolution operation, 
and the p phases of d( mT.) per sample period, T = pT, are 
described by 


‘d,=Sx,h,,-'u, i=0,---,p—-1, (46) 
, | 


where the ‘h, are i independent “subchannels.” With minor 
algebra, one can reduce the least-squares identification 
procedure to p subprocedures that can all be solved | 
separately. The p solutions can be interspersed to obtain 
Wi,,= W,,,, where N = M/p (we assume that p divides M 
or that M 1s picked slightly larger so that it does). Then, only 
one N x N matrix need be inverted (it is the same for all 
subchannels), rather than one M x M matrix, a considerable 
computational and storage saving. This matnx is the 
autocorrelation matrix of the input data alluded to in an 
earlier footnote (§). However, much greater savings are also 
available. 


e Use of fast algorithms 

The most efficient solution to the general least-squares 
identification problem appears in [2]. The DFT cannot be 
used in the general least-squares filtering problem because a 
Toeplitz structure must be imposed on (13) for its use. This 
solution requires about 


(2=+) IN + 4.5N? + pN° (47) 


multiplications, divisions, and additions in comparison to 
O(N’) for straightforward matrix inversion. [O(x) is a 
number that asymptotically rises no faster than in direct 
proportion to x.] The term (/N/p) + 4.5.’ is the fixed cost 
of the equivalent of inverting the matrix R,, ,, (fixed because 
it is the same for each subchannel); the remaining term 

pN’ + IN is the additional cost, at N° + IN/p per subchannel, 
for computing the equivalent of the product joe Pyyps* 
'W yp for each of the p subchannels. The storage 
requirements are about 6N + 2/ locations for the algonthm 
in (2]. The cost reductions accrue to the shifted nature of 
X,., with respect to X,.,_,, or equivalently, that R,.,, can be 
rewritten as a product of Toeplitz matrices, 


Ry ip = X sypsXvups > (48) 
where X,.,, iS defined in (26b). For more details, see [2]. 
e Choice of the input sequence 


Further computational and storage reductions are possible if 
the length-/ sequence x, , 1s chosen beforehand for all storage 
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=annels to be identified. A currently popular choice is a 63- 
“1 pseudorandom sequence. When the input data sequence 
:s Known beforehand, many of the quantities in the Fast 
‘53 TF) algorithm of [2] can be precomputed and stored 
once. reducing computation to 


me ae IN (49) 


inultiphications and additions (no divides) and storage 
(random access) to about 


Eo yatta (50) 
locations. Neither these counts nor the counts in (47) and 
(49) can be matched by the DFT or other methods of 
comparable estimate quality for reasonable N (20 or less). 
Asymptotically, because of the N log, WN computation in FFT 
implementations of the DFT, these FFT methods may have 
an advantage in terms of computational requirements, but V 
Is never chosen that large in practice. 


© fxperimental results 
lo demonstrate the robustness of the new least-squares 
identification method, several channel pulse shapes are 


IBM J RES DEVELOP VOL. 30 NO. 3) MAY 1986 ‘BM SJ 433 
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plotted in Figures 3(a) and 3(b), while the corresponding 
steps are plotted in Figures 4(a) and 4(b). These responses 
were obtained using the new procedure for a 63-bit 
pseudorandom sequence on digitized measurements of a 
thin-film disk/thin-film head channel [Figures 3(a) and 4(a)], 
and on a particulate disk/thin-film head channel [Figures 
3(b) and 4(b)]. The measurements were taken at several 
different diameters on each device. The diameters for Figures 
3(a) and 4(a) were 105, 120, 135. 150, and 165 mm, while 
those for 3(b) and 4(b) were 103, 136, and 172 mm. Figures 
5(a) and 5(b) show the pulse response and its spectrum, 
respectively, for an optical storage device. In Figures 6(a) 
and 6(b), we have plotted pulse and step responses for a 
magnetoresistive head in a magnetic-tape system; this time a 
62-bit pattern corresponding to NRZI coding of two cycles 
of a 31-bit pseudorandom data pattern was used [10]. In 
Figure 7, the delay for the magnetoresistive head 1s plotted 
to illustrate the ability of the new least-squares identification 
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Appendix A: Arbitrary sampling rates 
In this appendix, the sampling intervs! 7, is permitted to 
take the values 


q : 
as (Al) 


where g and p are relatively pnme positive integers such that 
q <p. Any arbitrary ratio of sampling to data rates can be as 
closely approximated as desired by the relation in (A1), as 
long as it is known, which implies some synchronization 


a between digitizer and write ciock. We also define a smaller 
3 time interval 7 by 
E b Tl, YT 
hes "Bal “ 
0.10 or 
qpt = pT, = qT. (A3) 
0.05 


The samples at rate 7, can be organized into successive 
3 disjoint sets of p members and of duration p7, = pqr. Then 
any sampling instant 7_, can be rewritten as 


mT, = npgr + iqr = (np + I)T, 


—0.10 1=0,---,p-l. (A4) 
0 40 80 120 160 
The equivalent of (37) becomes 


Time iterations 


d[npqr + iqr] = & hA(npqr — kpr + igr)x, 
; k 


+ u(npqr + igr). (A5) 


Note that, if p and gq are relatively prime. as was assumed 
earlier. then /: will be specified at intervals of 7 in (A5), or 
equivalently at all time instants that are integer multiples of 
7. At sample / within each group of p samples, only the 
procedure to capture that quantity as well. The dc level was —- Values (pr + iqr), where & is any integer, contribute to 
removed from the desired signal for the optical device to 
facilitate inspection of the plots; the true optical channel is a 
baseband channel. The plots in Figures 3, 4, and 5 
demonstrate the robust utility of the least-squares procedure. 400 


500 


4. Conclusions 

In this paper, we have introduced a new least-squares 
storage-channel identification procedure. We have analyzed 
the procedure thoroughly and demonstrated via expernment 
its utility and its improvements over existing methods. 
Several methods for reducing the implementational cost of 
the procedure were also discussed. The procedure can 
become a uniform standard for identifying and comparing 
the channel charactenistics of various storage media. | 0 0.2 0.4 U.6 0.8 1.0 
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a npqr + igr). Thus, d has p phases per group of p samples 
that can be independently modeled as 


d,, = p? X Angee ~ Up §= O05 D- |], (A6) 
where, again, 

‘h, = h(nT + iT,) (A7) 
and 

‘d,=dnT,+ iT,); ‘u, = u(nT + iT,) (A8) 


fori = 0. ---,p— 1. Each of the subchannels can be 
identified independently and the resultant responses overlaid 
(with delays of 7 with respect to one another). The overall 
response can then be used directly or decimated to pr (the 
data rate), gr (the sampling rate), or any other integer divisor 
of the rate 1/7. An important point to note is that there 1s a 
loss in resolution of a factor of approximately g for any fixed 
data length / with respect to the case where 7, = 7/p. This 
iast fact makes the alternative of resampling the data or 
phase-locking the ADC used to acquire the data (set g = 1) 
very attractive from a performance viewpoint. 


Appendix B: On-line efficiency 
It is possible to implement the least-squares storage-channel 
identification method in a sample-recursive manner. The 
procedure becomes a special case of the one considered 
previously by this author for echo cancelers in data 
transmission in [8]. The storage identification procedure 
could be performed on line, for example, to initialize, and 
possibly update (see [15, 16], a Maximum-Likelihood 
Sequence Detection Circuit. 

A brief summary of the procedure 1s, where k is the 
recursive time index, 


, r P 
Wan = Wagan + Ome + Cae » (BI) 
and where 
Ei = dk) = ere. Gre . (B2) 


C,,, 1s an M x | function of the input (presumably known 
or “training”) data sequence and is given by 


k 

ne aa) Me (B3) 

m=0 
and is presumably precomputed and stored prior to use. For 
more details on this procedure, and for an efficient recursive 
computation of C,,, when there is no prespecified training 
sequence, see [2. 8. 17, 19]. A final note 1s that. if the signal 
written just prior to the start of training 1s an erasure, then 
the prewindowed exact-initialization method of [8. 17] 
applies. rendering extremely low computational 
requirements: (B!) and (B3) simplify dramatically in that 


CdSe. 


IBM J. RES. DEVELOP. VOL 30 NO. 3) MAY 1986 ‘BM SJ od 


Appendix C: Methods for nonlinear identification 
The study ct nonlinear identification of a data channel is an 
entire subject area in itself. For instance, one can refer to 
[20] and [21] for methods based on simplification of 
Volterra senes under the constraints of a binary input. Here, 
a simple method suffices to verify the presence/absence of 
appreciable nonlinearities and to roughly quantify their 
magnitudes relative to the linear component of the channel 
response. 


SNR measurement 

Estimation of the SNR was discussed earlier. The minimized 
sum of squared errors, é,,,, contains a component due to 
modeling error. If 7 is sufficiently large. most of this 
modeling error is due to nonlineanties. The size of the SNR 
is indicative of the level of nonlinearities. Generally 
speaking. SNRs well below those expected can be indicative 
of large modeling errors due to nonlinearities. Thus. one can 
use the size of the SNR as an indicator of nonlinearities, 
given that he has some prior experience with the particular 
media and head and knows what to expect in terms of a 
nominal SNR value. This type of procedure requires a verv 
accurate phase-lock to the underlying data rate to ensure 
that nonlineanties are not artificially inserted by sampling- 
phase errors in the measurement process. 
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APPLICATIONS OF A PEAK DETECTION CHANNEL MODEL 


P.H. Siegel 


Abstract - A computer model of a peak detecting magnetic recording 
channel has been implemented and used for channel design and perform- 
ance evaluation. The model predicts raw error rate, ontrack and off- 
track, as a function of linear density, run-length-limited (RLL) modu- 
lation code, write precompensation rules, and tapped-delay-line (TDL) 
equalizer. It assumes noise additivity and validity of linear superposi- 
tion, and it bases calculations on a measured disk/electronics noise 
spectrum and digitized isolated transition readback signals from the data 
track and adjacent tracks. Details of the model are described, and 
illustrative applications to RLL (d,k) code selection and pulse slimming 
equalizer design for a specific channel are discussed. 


INTRODUCTION 


There are a number of signal processing options available which 
have the potential to increase areal density and reliability of peak de- 
tecting magnetic recording channels. Among these are modulation 
coding, write precompensation, and pulse slimming equalization. Assum- 
ing additivity of disk/electronics noise and adjacent track noise, and 
validity of linear superposition in the readback process, we have devel- 
oped a model of a peak detection channel which predicts raw error rate 
as a function of linear density and specified signal processing. The basic 
methodology employed is similar to that suggested by Katz and Camp- 
bell [1]. 

-Novel features of the model, in addition to the implementation 


exible write precompensation rules, and the incorporation of 
a TDL equalizer for general read equalization capability. 

We discuss below some of the technical aspects of the model. 
We then address two applications to a specific disk channel: a compari- 
son of RLL code performance, and the selection of a minimum noise 
pulse slimming equalizer. 


INTERSYMBOL INTERFERENCE AND CODE PATTERNS 


Intersymbol interference (ISI) affects the peak position and peak 
amplitude of the pulse resulting from a given transition. We compute an 
odd ISI length L, where L is the number of bits needed to account for 
pulse interactions. We calculate a cubic spline fit of a digitized read- 
back pulse from the data track, as shown in Fig. 1. Then, using linear 
superposition, we simulate the readback signal corresponding to each 
(d.k) pattern of length L having a central transition. The differentiated 
signal is also calculated with the spline coefficients. The position of the 
central peak is located by use of a Newton-Raphson iterative search for 
a zero-crossing in the differentiated waveform, and the central peak 
amplitude is then found. The model next computes the values of the 
differentiated waveform at the edges of the detection window corre- 
sponding to the central transition. Two types of clocking are consid- 
ered: an absolute clock and a mean-centered clock. The window for the 
mean-centered clock is centered around the average peak position 
described below, and represents an approximation to the window found 
in a channel with a PLL (phase locked loop). The waveform derivatives 
at the detection window edges are required for the bit shift error rate 
prediction. The average peak position is found by weighting the calcu- 
lated peak positions for all patterns according to the Shannon pattern 
probabilities, and summing. The pattern probabilities indicate the 
frequency of occurrence of each pattern in encoded random data for an 
ideal (d.k) code. Since run-lengths are uncorrelated in an ideal code, 
the pattern probabilities are found by taking suitable products of run- 
length probabilities which we compute using techniques from informa- 
tion theory. 


Manuscript received June 16, 1982. 
The author is with IBM Research aiberatons San Jose, California 
95193, U.S.A. 


readback signal. 


Write precompensation rules can be specified in order to reduce 
the effects of intersymbol interference. The rules are pattern dependent 
adjustments of the recorded transition positions: a transition is advanced 


or delayed by a specified amount according to the code pattern context 
in which it occurs. 


20 z 
E 30 i 
w = 
§ 20 2 
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E & 
= 10 3 
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Time (ns) Frequency (MHz) 
Fig. 1. Digitized isolated transition Fig. 2. Digitized disk/electronics 


noise power spectrum. 


NOISE STATISTICS 


The error rate calculation also requires a probability density 
function for noise and differentiated noise. For the disk/electronics 
noise, we digitize a noise power spectrum measured on a spectrum 
analyzer, as shown in Fig. 2. Normal probability plots of noise sample 
values measured from a dc-erased disk indicate that a Gaussian distribu- 
tion fits the data out to at least three standard deviations. We take a 
Gaussian distribution for the disk/electronics noise, with mean zero and 
variance given by the numerical integral of the measured spectrum. A 
Gaussian distribution for the noise leads to simplifications in dealing 
with the differentiated noise as well. The derivative, n’, of 2 Gaussian 
noise process n is again Gaussian [4], and the power spectrum T(f) of 
the differentiated noise is related to that of the original noise spectrum 
S(f) by the expression: 


TU) = (20f)” SW. (1) 


From the digitized S(f), we then compute the variance of the differenti- 
ated noise as the numerical integral of T(f). We model the distribution 
of n’ as Gaussian with mean zero and with this variance. 


Near track 


Far track Convolved 
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Fig. 3. Digitized adjacent track 
readback pulses, 4p offtrack. 


Fig. 4. Sample histograms from 
simulated adjacent track waveforms. 


For adjacent track interference (cross-talk), we assume that 
side-reading of the nearest adjacent track on each side of the data track 
dominates the cross-talk signal. The readback signal from an_ isolated 
transition written on the adjacent track is digitized for the head position 
of interest. Figure 3 shows the readback pulses from the near and far 
adjacent tracks when the head is 4u offtrack. Using a cubic spline fit 
and linear superposition, the readback waveform from several hundred 
bits of a pseudorandom (d,k) coded sequence is simulated and sampled 
up to 20 times per clock period. A histogram is made from the sample 
set as an estimate of the distribution density of samples from the adja- 
cent track. If we assume no correlation between signals from different 
tracks, the total cross-talk density is estimated by taking the discrete 
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convolution of the histogram densities from the two tracks. Histograms 
from a 4p head offset at a linear density of 18 kbpi with the (2,7) code 


are shown in Fig. 4, along with their convolution. The analogous calcu- 
lation is then carried out for the samples of the waveform derivative. 
The distribution for combined disk/electronics and cross-talk noise can 
then be computed by discrete convolution. 


ERROR RATE CALCULATION 


Given that the signal has a pulse peak in the detection window 
W, the probability of noise-induced bit shift error. is the probability that 
the differentiated signal plus differentiated noise waveform will fail to 
have a zero crossing within W: 

Pri s'(th + n'(t) < Oors'(th + n'(t)>0, fortinW). (2) 


Solving for this level-crossing probability exactly is a difficult mathemat- 
ical problem, even when n is a Gaussian process. We use instead a 
convenient approximation, suggested by A. Milewski, which is a reason- 
ably tight upper bound in the case of bandlimited noise. The probability 
is bounded above by the sum of the probabilities of the two events, for 
which simple upper bounds exist. For the first event, let ¢; be the time 
where the signal derivative is positive and of largest magnitude. Then, 


Pros'(t) + n'(t) <0, fortin W) § Pr(n'(t,) < -s'(t,)). (3) 


Similarly, if ¢, is the time where the signal derivative is negative and of 
largest absolute value, 


Pr( s(t) + n'(t) > 0, fortin W) < Pr(n'(t,) > -s'(tg)). (4) 


In practice, ¢, and ¢, have been found to lie at the detection window 
edges for the channels and densities studied. So, the approximations are 
evaluated at those points, using the waveform derivatives at the window 
edges and the noise distributions described above. We note that this 
upper bound has proven to be tighter than the approximation suggested 
in [1] which extrapolates the waveform derivatives at the window edges 
from the zero-crossing /, along a line of slope s"(tg) pee Fig.5. 

The probability of missing bit error depends on the clip level C, 
which represents the minimum amplitude necessary to detect a peak in 
the readback signal. The probability of interest is 


Pr( s(t) + n(t)| < C, fortinW). (5) 


This represents a level-crossing probability which we approximate with 
the simple upper bound 


Pr( |s(tg9) + Ut) 1 < C). (6) 


This probability is evaluated with the computed signal peak amplitude 
s(fg) and the noise distributions. See Fig. 6. 

These error probability bounds may be used for approximate 
worst case pattern analysis. A weighted average using Shannon pattern 
probabilities provides an estimate of the overall error rate for encoded 
random data. 


P. < Prin’(t,;)<-A,) ° 
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Fig. 6. Missing bit error rate 
approximation. 


Fig. 5. Peak shift error rate 
approximation. 


APPLICATIONS 


We now discuss two applications of the mode! to a disk channel 
with a thin film head and particulate medium. Measurements were made 
at the inner diameter. Track pitch was 30n. 
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RLL code comparison 


We predicted the performance of (1,8) and (2,7) codes at three 
head positions - ontrack, 2y offtrack and 4 offtrack. The clip level 
was assumed to be 40% of the base-to-peak amplitude of the data track 
pulse. No write precompensation or pulse slimming equalization was 
used. The resulting error rate/linear density tradeoff curves are shown 
in Figs. 7 and 8. For the (1,8) code, peak shift errors dominated at 
densities less than 20 kbpi, while missing bits were the major error 
mechanism at higher densities. For the (2,7) code, however, peak shift 
errors were the primary determinant of error rate at all densities consid- 
ered. The model indicates that at densities less than 20 kbpi, the (1,8) 
code has lower average error rates than the (2,7) code. At higher 
densities, the loss of signal amplitude degrades the (1,8) performance. 
In the range of error rates from 1E-12 to 1E-8, the (1,8) code provides 
a density advantage of slightly more than 5%. This result is consistent 
with the conclusions reached by Fisher and Newman in [5]. 

Table I shows a list of worst case patterns with L = 15 for the 
density 18 kbpi, as calculated by the model, along with peak shift, peak 
amplitude, and probability of error for ontrack operation. In general, 
the worst case patterns highlight features of the digitized pulses and can 
be used to assess the impact of peculiarities of pulse shape on error rate. 
Here, patterns with a minimum length run followed by a long run clearly 
affect the performance most severely, reflecting the pulse asymmetry. 


Log Error Probability 


14 16 i 20 22 24 
- Linear Density (kbp!) 


Fig. 7. Simulated performance of (1,8) code. 
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Fig. 8. Simulated performance of (2,7) code. 
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Peakshift Relative Probability Pattern 
(ns) Amplitude of Error 
tal 85 4.13E-5 100001010000101 
7.0 84 3.14E-5 010001010000101 
7.1 86 2.94E-5 100001010000100 
TZ 83 2.59E-5. 000101010000101 
7.0 86 2.23E-5 010001010000100 
(1,8) code, 25.1 ns window 
5.6 .90 1.73E-3 010010010000001 
eS 89 8.83E-4 000010010000001 
5.5 91 7.40E-4 100010010000001 
5.6 .92 5.77E-4 010010010000010 
a2 87 4.22E-4 010010010000000 


(2,7) code, 18.9 ns window 


TABLE I: Worst case patterns at 18 kbpi. 


Pulse slimming equalizer evaluation 


Barbosa [6] has reported on a design method for minimum noise 
pulse slimming equalizers. For a given channel and linear density, he 
constructs a One-parameter family of TDL equalizers, each of which 
maximizes the degree of slimming subject to a noise penalty constraint. 
At densities from-14 to 24 kbpi, we used the model to select the noise 
penalty for which the corresponding equalizer gives the smallest average 
ontrack error rate. The (2,7) code was used, and no cross-talk was 
considered. The ontrack and offtrack performance of the selected 
equalizer was then calculated, with cross-talk included. The results for 
the equalized channel with (2,7) code are shown in Fig. 9. 

The conclusion based on the ontrack performance is that these 
equalizers can increase linear density between 10% and 20% in the 
range of ontrack error rates from 1E-12 to 1E-8. The equalized channel 
is not sensitive to small offtrack excursions, but the offtrack perform- 
ance deteriorates as offtrack distance increases from 2p to 4p because 
of the enhancement of the cross-talk signal by the equalizer. 

The worst case patterns were found to reflect the positions of 
the sidelobes of the equalized pulse. For example, at 20 kbpi, with a 
detection window of 17.05 ns, and with the equalized pulse shown in 
Fig. 10, the worst case patterns had runs of 4 zeros preceding and 
following the central transition, thatis, 10000100001. 


= 
5 
4) 
i Equalized 
i 240 4u offtrack 
= ~~ Equalized 2u offtrack 
= 5 Equalized ontrack 
O 
ac 

—20 

14 16 18 20 22 24 


Linear Density (kbp!) 


Fig. 9. Simulated performance of TDL equalizer with (2,7) code. 
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Normalized Amplitude 
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Fig. 10. Comparison of unequalized 
pulse and 20 kbpi equalizer pulse. 


CONCLUSIONS 


We have described a computer model which predicts raw error 
rates for a peak detecting magnetic recording channel. Offtrack per- 
formance is predicted by inclusion of adjacent track interference effects. 
Calculations are based on measured channel characteristics: step respon- 
ses from the data and adjacent tracks, and a disk/electronics noise 
spectrum. The model also permits the evaluation of several signal 
processing options, individually and in combination: RLL code, write 
precompensation, and pulse slimming equalization. In addition to error 
rates, the model provides useful information about dominant error 
mechanisms and error sources both for worst case code patterns and for 
random coded data. Results of (d,k) code comparison and TDL equaliz- 
er evaluation for a specific disk channel were discussed. 
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Peak Shift Caused by Gaussian Noise in Digital 
Magnetic Recording 


Yasuhiro Tahara, Yoshimasa Miura and Yoshinori Ikeda 
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SUMDIARY 


The peak shift generated in digital magnetic 
recording processes is one of the most important 
obstacles to high-density recording. The princi- 
pal causes of peak shifts are waveform interfer- 
ence effects and noise. Of these causes, only the 
noise components have been subjected to empiri- 
cal treatment. In this paper we developed a prob- 
abilistic analysis of the peak shift due to noise. 
Gene rations of the peak shift are treated as prob- 
abilistic distributions and the corresponding dis- 
tribution functions and contribution to the phase 
margio are theoretically derived. The results 
show that when Gaussian noise is superimposed 
oo read-out signals from the head, generation of 
peak shifts due to the noise also exhibits a Gaus- 
sian distribution. With the variance of the dis- 
tribution as o¢, the maximum peak shift is 55 
~ 7 o and the loss of phase margin is ll~l4ic. 
The theory is applied tothe MFM recording sys- 
tem and the peak shifts of 2F, 1F and 1110] pat- 
terns due to the white noise are obtained. The 
ratio of these peak shifts takes an almost con- 
stact value of 1:1.3:1.2 In the region where the 
resolving power is 50 to 70%. It is found that 
the theoretical prediction and the experimental 
data agree very well for |110] patterns. 


1. Introduction 


Improvement of recording density is an im- 
portant problem in magnetic devices such as mag- 
gpetic disks and drums. One of difficulties en- 
countered in high-density recording systems is 
the generation of peak shifts. Information writ- 
ten on disks and drums as magnetic reversal pat- 
terns is read out at magnetic heads and regeao- 
erated by detecting the peaks of read-out wave- 
forms. When the recording density on disks and 
drums is increased, peak sbifts in the read-out 
waveforms become larger due to interdit waveform 
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interference, resulting in degradation of the SN 
ratio. If peak sti®s become excessively large, 
read-out errors occur and it is no longer possible 
to derive recorded informations from the read- 
ow waveforms. 


Principal causes of peak shifts are waveform 
interference and waveform jitter due to noise. 
A number of theoretical and experimental studies 
have been conducted on the waveform isterferenoce 
effects [1-3] and some efforts to reduce peak 
shifts have been tested using circwit technicues, 
such as wavefor=: equalization, that make ase of 
waveform characteristics [4, 5]. On the other 
band, only empirica! treatments have been done 
on peak shifts caused by noise and no specific 
Quantitative analysis has been conductec. In 
present-day magnetic recording devices, SN 
ratios are steadily decreasing because of the re- 
duction of read-out vaitages in the head as accom- 
panied with increased recording densities ard be- 
cause of the increase in noise Sandwicth due to 
higher recording and regenerating frequencies. 
As aresult, the effect of the noise has increased 
and methods must be developed for quantitative 
analysis of peak shifts caused by the noise. Fur- 
thermore, a design procedure for recording and 
regenerating syStems is needed by whicn the over- 
all peak shifts caused by both waveform interer- 
ence and noise effects may be minimized. 


One of the characteristics of the peak shit 
caused by the noise is its rancomness. This is 
because the noise generation is also random. 
Heace, for quantitative treatment of the peak 
shifts caused by nolse it is necesSary to introduce 
a probabilistic approach. Although Melincsoa [6). 
and Kobayashi [7] analyzed noise in NRZ record- 
ing systems in a probabilistic maaner, toey bave 
not considered peak shifts at all. 


In this paper, a method previously proposed by 
the authors [§] is extended to the probabilistic 


analysis of peak shifts caused by noise, and the 
probability distribution of peak shifts and their 
contribution to the phase margin are quantitatively 
derived. The theory (s applied to the case of 
BIFM recording systems, and the amounts of 

peak shifts caused by the aoise are computed for 
several practical patterns. Finally, experimental 
results are compared with theoretical predictions. 


2. Peak Shifts Caused by Noise 
2.1 Noise 


In magnetic disks and drums, information writ- 
ten on the recording medium as magnetization re- 
versal patterns is read out at the magnetic head 
and immediately amplified by a preamplifier 
located near the head. The amplified signal ls 
then sent to a peak detector at the later stage of 
the ae Principal causes of noise are (see 
Fig. l: : 


1. Preamplifier noise 
2. Head impedance noise 
3. Medium noise 


The first of these causes arises from semi- 
conductor noise generated in the preamplifier and 
consists of thermal noise, shot noise and 1/f 
noise, At the 1 to 10 MHz used in magnetic disks 
and drums, 1/f noise Is negligible, anc hence 
thermal noise and shot noise are predominant. 
The spectrum of the latter two is almost Identical 
to that of white noise. The bead impecance noise 
is a kind of thermal noise caused by the composite 
impedance of the head and head termination circuit 
as seen from the preamplifier, The spectrum 
distribution is concentrated near the resonance 
frequency of the head [9]. The medium noise is 
caused by nonuniform dispersion of magnetic 
particles in the recording medium (10, 11). Bis 
read out with the information signal by the head. 


These kinds of noise are generated randomly 
and hence the noise distribution can be treated as 
Gaussian, 


Let us expand the coise voltage Vy(t) into Four- 
lerseries at-T stsT 


W(t)= Fee sings + d, cos %,t) (1) 
om 
where 
e,= 2zf,= =xn/T (2) 


and ap and bp are probability variables that inde- 
pendently obey Ganssian distributions, Distribu- 
tion functions of ag and bp are identical and their 
mean values are zero. 


(2 


The power spectral density of this noise (de- 
fined as the mean square noise voltage per unit 
bandwidth at frequency of f,) is given by 


5 (a + 4,2 ) 
Sf 
where Af =1/2T. 


(4) 


NGU,)= 


2.2 Peak shift caused by noise—sinusoidal 
waves 


Let us coosider what kind of peak shifs will be 
produced when the noise described above is super- 
posed on the regenerated signals, We first ex- 
arnine the most fundamental case, in which the 
read-out waveforms are described in terms of 
Sinusoidal waves. The total regenerated signal 
valtage V(t) is 


Vitv=sy, cos t+ ) (a,sin ays +d,cos 2) (5) 
Pe | 


where Vo and wo are the amplituce and angular 
frequency of the regenerated signal, respectively. 


We shall now obtain the shift of the peak, origi- 
nally located att = 0, caused by the noise. Fir” 
V(t) is differentiated and expanded around t = 0, 
assuming the amour of the peak shift is small, 
The result is 


d ] 2 
ome = ea 6 
mF V(t) 5 2 Vit J en 2 (6) 


From this equation the amount of the peak shift Is 


2 
s= er 
Uy Vo (7) 


The mean square value of oc of At Is given by 


2 
4 ee 
of = = w, a, 
( St y Gy‘ Ve" e n @ 
— Fw? a2 
Les ee 
4 T 3 
= 4 2 ya NA) Sf 
Ve qm | 
] p= 
=e Te v(w)co 8) 
Gy Vy 20 


where N(w) = N(fq). 
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Fig. 1. Noise spectrum in magnetic recording file. 


$s Vo (t) 


Fig. 2. Read circuit. 


Next, let us examine the distribution function 
of the peak shift At. Notice that ap in (7) has a 
Gaussian distribution. In general, when variables 
xi (i =1, ..., N) have independest Gaussian dis- 
tributions the variable x given by 


(9) 


A 
= y c, =, 
a* 1 


where the cj are constants, also has a Gaussian 
distribution and its variance is given by 


y 
2 92 
be? ¢ 


a= 


02 = 


(10) 


where o|@ is the variance of x; [12]. Hence, if 
@q bas a Gaussian distribution, so does At. The 
variance of At is of course given by (8). The dis- 
tribution function p(At) is now 


2.3 Peak shift caused by the noise—geseral 
case 


1 (3)? 
/2% 6 aaa 20? 


where c~ is given by (&). 


P(Ji)= (12) 


Let us now derive a method for computing peak 
Shifts in general cases. The read circuit of a 


- 


9 


digital masnetic recordirg device generaily con- 
sists of an amplifier, low-pass filter, differertia- 
tion circuit, and zero-croussing detector as stown 
in Fig. 2. Since operation of these circuits is not 
ideal, the effects of their frequency characteris- 
tics on the peak shift cannot be neglected, 


For simplicity, we represent the frequency 
characteristics of the entire read circuit by that 
of a low-pass filter. When the transfer function 
of the low-pass filter is F(s), the transfer func- 
tion H(s) of the read circuit, which cantains a dif- 
ferentiation circuit, is giveo by 

H(s)=seF(s) (12) 
as is seen from Fig. 2. If an RC approximate dif 
ferentiation circuit shown In Fig. 3 is used in plac 
of a true differentiation circuit, the tracster func: 
tion H(s) becomes 


s 


F(s) (13 


s 
Is )=——- F(s)= i 


w, 


S 


wy 


where wq =1/RC. Hence, if another first-order 
low-pass filter is inserted, the transfer function 
can be reduced to that of (12). 


Consider now the case in which a read-out Sig: 
nal f(t) is indicext at the circuit, Since the total 
input signal including the noise is 


Fig. 3. Approximate dif- 
ferentiation circuit. 


V(tj=f(t)+ y (¢,sinw,t+ bcos?) (14) 


the output Vpi{t) of the differentiation circuit is 
Signal a 
Vol t)=sle)+ Pas |FOe,) | 
X{ 4, cos («6+ 6.) = sb sin (o¢+ 9 )} (15) 


where h(t) is the output signal of the differentiation 
circuit for the read-out signal f(t) and is given by 


LCAC*) = F(s)- LU (2)) (16) 


using the transfer function H(s). 
of F(jwy). 


On is the phase 


Let the zero-crossing point of h(t) be to and 
the slope of h(t) near to be Go, lL.e., 


d 


Then the zero-crossing point of the output Vpit) 
is given by 


(Qt ps we | F(a) | 


X{@, cos (“sf +0,) - 5, sin(% 4 +0n)} (18) 
The first term to represents the peak shit caused 
by the waveform interference and the phase delay 
in the circuit, whereas the second term corre- 
sponds to that caused by the noise. Hence, the 
mean square value of the peak shift caused by the 
noise is 


ie 


7s ooh oF | F(Jw) | vw) dw (19) 


or, using (12), is given by 


0? = =a aes |e! N(w)dw (20) 


where the term 


sof Lue) | v(w)de 


represents the ncise power contained in the output 

of the differeoriation circuit. Therefore, the peak 
shift caused by the noise can be described in terms 
of the noise power in the output and the slope at the 
zero-crossing point. 


3. Peak Shifts of Various Patterns in ATM 
Recording Systems Due to Noise 


In conventional digital magnetic recordings de- 
vices, the recording and read-out of information 
are performed using PM, FM or MFM recording 
processes in which self-locking can be incorpor- 
ated. In most recent large-capacity recording de- 
vices, the MFM process is employed. In this sec- 
tion, the effect of waveform interference on the 
peak shift caused by noise is investigated. To this 
end, peak shifts caused by white noise are calcu- 
lated and compared for 2F, 1F and {110} patterns 
(Fig. 4). 


3.1 2F patern 


In the MFM recording processes, the 2F pat- 
terns have the highest magnetization reversal fre- 
Quency. Ordinarily, the recoruing density in the 
MFI process is such that the resolution is 60 to 
70%, within which range the read-out signals of the 
2F patterns are almost sinusoidal and the contri- 
bution of harmonics is negligible. Hence, the 
read-out signal waveform f(t) is 


ve (£)=F Ky ecos aot (. 


where VoF {s the peak-to-peak amplitude. Here 
WO (= 27fq) is the recording and read~out anguler 
frequency and is related to the bit cell time T via 
v= x/T (22) 

Let the power spectrum of the white noise be 
N(w)=n¢ (23 


We assume that the low-pass filter Is an Ideal with 
the characteristics 


1 isSsf 
F(jw)= | J (24) 
0 if>t. 
Then the peak shift is computed from (19). 
Y2, eg 
SS aw 2s 
a v3 xf? Voy 
If the noise power is written as Nrms- 
Nos= J, : 7 (26) 
then OOF 
te N 


Co = owe GQ 
3F 
V3 74,7 2F 
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Fig. 4. Read-out waveforms of various patterns in MFM recording, 
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Notice that the amount of the peak shit is pro- 
portional to {.3/2, as shown In (25). This is be- 
cause the higher frequency content of the noise 
contributes more to the peak shift. Hence, inthe 
read-out of magnetic recording systerss the im- 
portance of elimination of high-frequency noise by 
the low-pass filter is more than just for improve- 
ment of the SN ratio. 


3.2 1F patterns 


The 1F patterns have the lowest magnetization 
reversal frequency in MFM recording processes, 
As seen from Fig. 4 (b), the read-out signals 
contain the 3F harmonics. Letting the amplitude 
of the read-out signal be ViF and the ratio of the 
3F components be 8, we can write the read-out 
waveforms of 1F patterns as 


fel(th=z Ke{ (1A )eos 5 a, + cos 5% ! (28) 


where # is usually 0.1to0 0.2. The amount of 
peak shifts of 1lF patterns caused by the white 
noise can be obtained in a manner similar to the 
case of 2F patterns and is given by 


: 4f, \ ems 
= 3xrf2 (1483) My 


The peak shift of 1F patterns depends on ViF 
as wellas 8. This is because the amouct of shift 
is not a function of the amplitude but of the sharp- 
ness of the peak. Hence, the peak shift caused by 
the noise is less likely to occur for the case with 
larger resolving power and 8. 


(29) 


Next, the magnitudes of peak shifts of 1F and 
2F patterns will be compared. If the noises are 
identical in both cases, the ratio of o1F and cor 
is obtained from (27) and (29) as 


AY 2 AR 
Oop 1+8 7 (30) 
where R is the resolving power given by 
R= Vis/Vip (30) 
Let an isolated read-out waveform be 
@ 
ei)=A—_— (32) 


at + (2 


where A and a are constants for expressing wave- 
forms. If we assume that the read-out signal 
waveforms can be expressed interms of a super- 
position of isolated waveforms, & is given by 


i-vi-R (33) 
2 


(see the Appendix). Hence, (30) becomes 


sp 4R 
mee sieeecemeeeieatic (34) 
6.9 §-4/1-# 


Figure § shows o1F/o2F versus R. From the 
figure, it is clear that 1F patterns are approxi- 
mately 30% more likely to be affected by noise 
30°. 


3.3 {110} patterns 


In MFM recording systems the maximum peak 
shif caused by waveform interference occurs for 
"110110" patterns. The read-out waveforms of 
these {110} patterns can be represented by using 
(4/QF and (8/3) F components as 


Sirol! YEdV,(sinZ oy ¢- r sing 04! (35) 


where V110 is the amplitude of (4/3)F components 
and y is the amplitude ratio of (8/3)F to (4/3) F 

- Components. Peaks of these {110} patterns are 
located at 


(36) 


_ 37 “(ee 


‘=—ee 
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and the amount of peak shifts AT 110 caused by 
waveform interrerence is 


3 Bs 1V1+327r?)\ } 
AT =T |) 08 ( ae (37) 


When white noise is superimposed on these {110} 
patterns, the amount of peak shift caused by the 
noise is 


; 3V3f _ rms __ (38 
po AR fz VAT), ) 

2 
Pn ie alee es al aes iN tee (39) 


2 leV1+732 77 ) 


As In the case of 1F patterns, we consider the 
tratio of C110 too2r. If we assume that (32) 
represents the isalated read-out waveform, we 
obtain (see the Appendix) 


4/3 
= P=" I (40) 
As 
Vy v 37 (41) 


Hence, if the noise is identical in both cases, we 
have 


Siw 2 in (42) 
Le 4 A(7) 
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Fig. 5. Relative strengths of the peak shifts 
caused by noise as a function of resolution. 


Figure § shows o110/CoF versus R. From the 
figure, we observe that {110} patterns are almost 
20% more likely to be affected by noise. 


The electromagnetic conversion characteristics 
of magnetic heads are often expressed in terms af 
the resolving power and the output voltage of 2F 
patterns. Act&zal peak shifts caused by the noise ~~ 
are larger for 1F and{110} patterns than for 2F 
patterns. For resolving power of 50 to 70%, the 
amount of shifts exhibits a constant ratio of 1:1.3: 
21.2. Therefore, the magnitude of phase shifts 
due to noise can be predicted from the output volt- 
age, and hence from the SN ratio, of 2F parerns. 


4. Phase Margin ard Error Rate 


The demodulation process in usual digital mag- 
netic recording systems consists of the fallowing 
steps. First, using a phase-lock loop, clock sig- 
nals are generated from the data pulse train 
emerging from the peak detector. From these 
clock signals, window pulses are then created. 
Discrimination of 1 and 0 is done by detecting 
whether a particular data pulse is within the win- 
dow pulse, The phase margin is defined as the dif- 
ference between the maximum peak shift actually 
generated and the width of the window pulse, This 
quantity ls viewed as a figure of reliability for mag- 
netic recording devices. When the phase margin 
is sufficiently large, correct demodulation is pos- 
sible even if new peak shifts are generated due to 
small defects on the recording medium or to track- 
ing errors, as long as their magnitudes fall within 
the phase margix. 


As shown in Fig. 6, measurements of the phase 
margin are performed by shifting the window pulse 
with respect to the data pulse and by detecting the 
error rate. The phase margin is the width of the 


distribution 
of peak shift. 


window pulse (nor!) 


window pulse (shifted ) 


Error rate 


p (4t) 


47 


| 
| 


Window position 


Fig. 6. 


Principle of the measurement of phase 


margin. 


window pulse shift for which the error rate Is be- 
low a certain value. 


Let us derive the probability P(AT) at whicha 
Gata pulse creates a peak shif larger than AT 
due to noise, The result is 


oe 4T 
P(sT)= J p(4t)dst=erfe (=) (43) 


where eric(x) is the erzor function. Figure 7 
shows P(AT) versus AT/co. Siace the error rate 
in magnetic recording devices should be 1078 

~ 107+, we see that a maximum peak shift of 
5.5 ~ 7o is created. The phase margin is re- 
duced by 11 ~ l4o due to the noise. 


In general, when a peak shift of ATk (k =1, 


eee, N) already exists due to waveform interfer- — 


ence, the probability of generating peak shifts 
larger than AT is 
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N 
B(4T)= P w, P(4T- 47) (44) 
ml 


where wk is the ratio of pulses which cause the 
peak shift of ATk. 


5. Comparison with Experiment 


In the experiment two kinds of magnetic heads 
and disks, A and B, were used to record and read 
out and the variation of the error rate was mea- 
sured with respect to the location of the window 
pulse. The recording process was MFM and the 
recording and read-out frequency was 6.45 MHz. 
The cutoff frequency of the low-pass filter in the © 
read-out circuit was 11.7 \idz. The read outpur 
of the head (Vor), resolving power and noise were 
measured and are listed in Table 1. Althoush the 
bead output of A is larger, so is the noise in A. 
The SN ratio of B is better by about 1 dB. Hence, 


Table]. Recording and read-out characteristics 
of A and B 


ee B 
Head output L32 mVp-p | LOS mV, , 
Resolving power | 61% 67 & 
Noise ) 312 2Vrms | 22.72 Vrms 


47ico=— 
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Fig. 7. Probability of peak shift larger than 
AST caused by noise, 
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it is expected that the resolving power and SN 
Yatio are better and the peak shift smaller in B. 


Figure 8 shows measured results of the error 
rate with respect to the window pulse location 
when {110} patterns are recorded and read out. 
As expected, the peak shift of A is much larger 
and the loss of phase margin is greater. 


Values in Table 2 were calculated from the 

measured values in Table l. In Table 2, the peak 
shift caused by waveform interference and the peak 
jitter due to noise are listed. Computed results 
of the error rate are solid lines in Fig. 8. Their 
agreement with experimeastal data is excellent, 
In the present experiment, the loss of phase mar~ 
gin due to nolse was 30 to 40% of the window pulse 
width and we can see that the effect of noise on the 
peak shift is quite important. 


6. Conclusion 


The peak shif caused by noise was analyzed in 
@ probabilistic manner and several examples were 
computed. The distribution and magnitude of the 
peak shift and its contribution to the phase margin 
were studied. Although, to date, peak shifis have 
been analyzed using empirical methods, the new 
method in this paper is capable of predicting more 


Table 2. Computed values of peak shifts 


Waveform interfer. 47,,, 
Jitter by noise °%,;, 


6.0ns 
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Window pulse position (ns) 


Fig. 8. Experimental results of the error rate for 


the {110} pattern. 
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accurate values. In future designs of magnetic 
beads or recording and reading circuits, the total 
peak shif caused both by waveform interference 
and by noise must be taken into account. The 
present method is useful for the optimum design 
of such devices. 


One of the problems yet to be analyzed is 
waveform equalization by such circuits as pulse- 
narrowing networks. When the read-out wave- 
forms go through a waveform equalizer, the peak 
shift due to waveform interference may be re- 
duced, whereas that due to noise may increase. 
Waveform equalization Is useful only when the re- 
duction of the peak shift due to waveform inter- 
ference Is larger than the increase of the peak 
shift caused by noise, Since the SN ratio gradu- 
ally decreases as the recording density Is in- 
creased, the design of equalizers most be done 
with extreme caution. 


In the present paper the peak shiff was as- 
sumed small in order to simplify calculation was 
simplified. When the peak shift is extremely 
large, this simplification {s no longer valid and 
the distribution is degraded from a Gaussian 
form. However, in conventional devices, the SN 
ratio is larger than 20 dB and the simplification 
oo above is believed not to cause any prob- 

ms. 
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APPENDIX 


Let the isolated read-out waveform be 


e({)= As 


@ 
a? 27 
and the bit cell time be T. The read-out wave- 


forms of 2F patterns in MYFA1 recording systems 
are 


Bpl= Fo CarseCe ne) 


Fourier transforming the above, we obtain 


faglt)= 240, LT exp{-( 201) aye} ms { (21) w, ¢ } 
Rand 


where 
o=2/T 


Slnce costributions of the second- and higher- 
order terms are Small, they are neglected 


Sygl 8)= 2 da, exp (-@, 2) wos w, f 
In the case of 1F patterns, contributions of the 


third=- and higher-order terms are similarly ne- 
giected, 


fy) = Aw, {exp (-4 Wo <) cost w,¢ 


+ exp (-$,«) cos 3 wot ! 


k= 240,0= exp{-tu, e) 


Then, V2¥F and ViF become 


Viy=2ke 
Vip=k (are?) 


From these equations, we obtain 


Vy 1+a? 


a? 1-/1-R? 


1+a? 2 


f= 


The read-out signal waveforms of {110} patterns 
are 


reatn | 


ry 
R 


From the above, we have 
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Configuration of Codes 


Error oe Modulation Modulation Error User 
Correction Channel Correction 
Code Code Data 
Code Code 
Encoders Decoders 
Modulation Code: . Matches recorded signal characteristics to channel 


bandwidth, detection method, read/write ‘electronics, 
timing and tracking servo requirements 7 


Error Correction Code: Detects and corrects data detection errots 
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Digital Magnetic Recording Channel 


Magnetic 
Track: 


Data: . O . 1T . 1. OO . NRZI 


Write 
Current: 


Readback _-2% ~ 
Voltage: | ~~ 


Detected 
Data: : O.-, 1 : 1, O 


Causes of bit detection errors 
e Random noise 
e Intersymbol Interference Bit pattern 
e Loss of clock synchronization related 
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Digital Optical Recording Channel 
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Run-Length-Limited (RLL) Codes 


e (d,k) constrained codes 
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state diagram representation of (1,3) constraints 


Code rate = (user bits/code bits) 
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Code Constraints and Capacities 


Unconstrained 


input 


|—m-| 
000000 
000001 
000010 


111111 


Shannon capacity = C < 1 


Constrained 


output 
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001001001 
010001001 
010010010 


100000001 


Necessary condition for block code: 


Sufficient condition for block code if n large: 
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Example: (d,k) constraints 


C = logo A, where A is largest real root 
of p(x) = xkttexk-d 
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Magnetic Disk 
1/2 (0,1) FM, Double 
Frequency, 


Manchester, 
Biphase 


1/2 (1,3) Modified FM (MFM) 
Delay Modulation, 
Miller 


2/3 (1,7) (Jacoby) 


1/2 (2,7)  (Franaszek) 


1/2 (2,11) 3PM 
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Examples of RLL Codes 
Magnetic Tape Optical Disk 
1/2 (0,1) Phase Encoding 


4/5 (0,2) Group Code 
Recording (GCR) 


8/9 (0,3) 


1/2 (1,3) Zero Modulation 1/2 (1,3) Delay Modulation 
(ZM) 


Miller2 2/3 (1,7) | (Horiguchi-Morita) 


8/17 (2,10) EFM 


Comparison of RLL (d,k) Constraints 
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(2,7) vs MFM (1,3) 


Fixed User Bit Density 
MFM a UW «& ODO =~ J 
(2,7) » tt. OO. OO. 1 


(2,7) reduces intersymbol interference 


Fixed Minimum Transition Spacing Twin 
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(2,7) 4 


(2,/) increases density ratio by 50% 
with decreased detection window 
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Tin 


Tmin = 3 


Density Ratio | 
2x 1/2 = 1 bit/T pip 


Density ‘Ratio 
3 xX 1/2 = 1.5 bit/T pip 


Peak Detection Channel Model 


Peak shift error rate calculation 
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(2,7) or (1,7) 
with 
enhanced channel 
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Linear Density (KBPI) 


Factor of 2.5 in linear density attributable to 


RLL code progress 


Power Density (db) 
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RLL Code Spectra or 
Fixed Data Rate CBA et wi 
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Spectral Null Codes 


Null at DC (f=0) Dc FREE 


— AC-coupled read/write electronics, 
rotary head (Het Scact) 


Power 
— Low frequency noise suppression Density 
— Tracking and focus servo on 
pregrooved disk 
I re) ice O 5 fo 
Ay PRRTCUAL FLEQVENCIES 
Nulls at f40 
— Embedded tracking servo 
(magnetic disk) 
ks Power 
— Embedded timing servo Density 


(optical disk) 


0 4,712 fo/2 ( 


DC Null Codes 


e Running Digital Sum (RDS): “accumulated charge” 


Bits .0O .0O .1 .0 .90 . 0 


Write | 
Signal 


Levels 
{aj} 


N 
2 a | <c, for all N = 1 
all {a} 
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e State diagram for bounded RDS signals 
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Maximum rate: logs 2 cos 6 . 5) 


e RLL combined with RDS: (d,k;c) constraint 
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Codes with Null at f#0 


e Generalized RDS at f=kfo/n 
N 

RDS,(N) a aie lemki/n = DFT (ao er an) 
i=0 

Spectral null at f<=>|RDS,(N)| < c, for N20 


e State diagrams: examples 


aa aaa 
LAAANW 


(5/2 


f=DC 

and 

f=fo/2 
(Interleaved DC) 
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Techniques 


Block codes 


Sequence state codes 
(fixed and variable length) 


Look ahead codes 


Sliding block codes 
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Code Construction 


Examples 


4/5 (0,2) 
8/9 (0,3) 
1/2 (0,1;1) 


1/2 (2,11) 
2/3 (1,7) 
1/2 (1,3;3) 
8/17 (2,10;c) 


1/2 (2,7) 
2/3 (1,7) 


GCR 
FM 


MFM 

(Franaszek) 
(Horiguchi-Morita) 
Miller? : 


3PM 
(Jacoby) 
ZM 

EFM 


(Adler-Coppersmith-Hassner) 


Practical Code Implementation 


Encoder Finite State Machine 


Modulation 
encoder 


(states) 


Features 


e High rate m/n 
e Low complexity 


Decoder Sliding-Block Decoder 
cere eok 
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decoder window (S\¢r REG} 


Modulation 


decoder PR. Oh Ole Ton 


Features 


e Limited error propagation 
e Low complexity 
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Block Coding 


FM (Frequency Modulation) 
1/2 (0,1;1) 


Encoder: Insert redundant code bit “1” between 
consecutive data bits 
(clock synchronization and dc-balancing) 


Decoder: Drop redundant code bits 


Example: 
Data: .1.1,. 0. 0. 
FM Coded: o11.11.10.10. 


Signal | LI lt 
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Sequence State Coding 
(fixed length) 


MFM (Modified Frequency Modulation) 


1/2(1,3) 

Encoder: 

A. After data bit = “O”: B. After data bit = “1”: 
0 ~~ 10 0 — 00 
1 57.04 ae Oe 


Data 
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Decoder: Drop redundant bits 


Example: 
Data: . 1. 1. 0. 0. 
MFM Coded: .01,01,00.10 
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MFM 


(Modified Frequency Modulation) 


1/2(1,3) 
Encoder: 
A. After data bit = “O”: 


0 ~~ 10 0 — 00 
bees OH 1 — 01 
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Decoder: Drop redundant bits 


Example: 
Data: . 1. 1. 0. 0. 
MFM Coded: .01.01.00.10 


B. After data bit = “1”: 
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| MFM 1/2(1,3) 


a | 
e 91% efficient (0.5/0.5515 =~ 0.91) | 
| 


e Finite-state fixed-length encoder (1 bit —2 bits) 


Encoded 
Data: Data | 


Pe baal 


os 


e Sliding block decoder 


State A: Previous input = “OQ” 
State B: Previous input = “1” 
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e Error propagation <!1 data bit | 


ne 


Data 


Clock 
Phase | 
= Odd 


State 
Indicator: 


Data. 
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state diagram G for (1,3) 


Construction of 1/2 (1,3) Code 


01 00. 


G 2 
1/01 0/00. | 
gO a 
0/10. 
1/01 
0/10 


Encoder graph 
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State-transition matrix T 
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1/01 
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Simplified (MFM) 
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Sequence State Coding 
(variable length) 


e 1/2(2,7) (Franaszek) 


e Graph representation of (2,7) strings 


e Finite-state encoder based on graph states 
— Fixed length code impractical (34 bit codewords) 


— Variable length block code 


0100 /C 
1000 /D 


000100 /C 
100100 /C 
001000 /D-> 
00100100/C 
00001000/D 


e Practical fixed length encoder obtained 
by introducing new states 


e Sliding block decoder with 
error propagation =< 4 data bits 
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(2,7) Code Implementation 


e Encoder logic circuit e Decoder logic circuit 


nein PLL Rae 
y, 


Indicator 
Encoded Encoded 
Data Data 


Om 


State 
Indicator | A 


Look Ahead Coding 


e 2/3 (1,7) (Jacoby) 


— Basic encoding table 


e Potential (1,7) violations: OlOloe! 
00.00 — 101.101 


— Violation substitution table 


rN 


101-0 00 
100:0 0 0 
001°000 
010-0 00 
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(1,7) Code Implementation 


- @ 98% efficient 


e Finite-state fixed-length encoder (2 bits — 3 bits) 


100/B | 010/V | 010/B 
100/C | 001/C | 010/C 
100/D | 001/D | 010/D 


State A: Previous input = “00” (no violation) 
State B: Previous input = “01” (no violation) 
State C: Previous input = “10” 
State D: Previous input = “11” 


State V: Previous input caused “Violation” pattern 


e Sliding block decoder with 
error propagation < 5 data bits 
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EFM 


e Eight to Fourteen Modulation 
e 8/17 (2,10;c) 
e Look-ahead encoder 


Data Bits ol itsfojofojoty of fofofotsfsyt, 
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RDS 


Decision Point 


e Block decoder (after discarding merge bits) 


e Error propagation < 1 byte 


Siding Block Code Algerithm 
{(Acdier, Coppersmith, Hassner, Niarcus) 


General code construction procedure for finite 
memory channels, e.g. (d,k) 


Produces code at any rate m/n < C 
7 Finite-state encoder 

¢ Sliding block decoder 

4 Limited error propagation 

Based on results in symbolic dynamics 


“Automatic” code construction possible 
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Sliding Block Code Algorithm 


e Generates new graph representation of (D,K) constraints - 


e Finite-state encoder based on new graph states 


e Example: 1/2(2,7) 
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State-Splitting 
21 states 
Two 2-bit codewords/state 
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State-Amalgamation 


7 states — 
Two 2 bit codewords/state 
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Encoder/Decoder Logic. 
for (2,7) code 
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Encoder logic circuit . \Decoder logic circuit 
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Constrained Codes for Partial Response Channels 


(0,G/I) constraints for (1 — D?) 


e 0 = minimum run of 0’s 
e G = maximum run of 0’s in channel output (Global) 


© | = maximum run of 0’s in even/odd substrings 
(Interleaved) kK Consrbinit ov Even PO C90 Subsreines 


Why? 


© 0 = no restriction on intersymbol interference 


e G = timing/gain control information 


@ | = limit iength of minimum distance error evenis 

Applications 

e 8/9 (0,4/4) and (0,3/6) block codes 

e 8/9 (0,3/5) sequence-state 
code 


© 8/9 (0,3/3) sliding-block code 
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MATCHED SPECTRAL NULL CODES 
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Improved codes for low SNR optical/magnetic’channels 


e Code design principle based on new general theorem 


null frequencies in 
code power spectrum 
au” que” eee 


EMM 
(rate 2/3) 


PR1:3 dB 
PR2: 4 dB 


Interleaved — 


Biphase 
(rate 1/2) 


PR4: 4.8 dB 
EPR4: 4.8 dB 


Significant coding gain results if: 


_ Jnull frequencies in t 
~ [channel frequency response 
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Outline 


Introduction and Overview 
Information Theory 
— State diagrams 


— Shannon capacity 


Statistical properties 

— Power spectrum calculations 

Design of Run-Length-Limited Codes 

— Code construction techniques 

— Applications 

Design of Spectral Null Codes 

— Characterization of spectral null constraints 
— Applications 


Trellis Codes for Partial Response Channels 
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Loops al lowed 
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Coding PRBLEM : Convert arbitrary (unconstrained) deter 
an shit nathan efficiently itp Sequences 
a“ -_ constTained by DNC. 

Defu: Th mh power , s™. of o DNOE LS is the 
system obTaimed by blocking SRopuences of S 
nto (non-overlapping ) Grou ps of ™ Symbols 


X eet KX Kye Ky. Xn croX 


x, xX, xX 


Remark ; TE Ss described by FLTD GC, Then SM 5 
deser bed by GC”, whi has Bi shes as G, 
and an aden ctate 4 tb tile} for onery 
poth of length me Aew Lt 4 in G, well corresponding 
m-block fabe- 
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Defw: A rate M:N Code from Ue, 4S is A pair of 
wa ps ‘ 
Ender. E: UL ih [Ee |— 
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whore Do bE = Tr, the pal Function 
it. D(EM)=x- 
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Finite- State -Mathine Encoder 


Def, A finste-etate- machine (FSM) is a £ynchronous System 
Consisting of : 1) a finite wiput alphabet X = S0%,,..% | 
2) & Finite output alphebet V = [fy 6x3 
3) a finite ctak set Z = £6,...,07, 
u) A pair of Cura cteriting bunctions 
Ty + ott function 
anh Fe: stete function 
aver _ ~ 
ouveus Ye = Ty (Xe Ze) 
Nexr sTh% Atay = f, (Xe, Z) fr. Came t =1, 2, us 


Graphical representation 
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S Tiding Block Decoder 


Def: D: he Us" iS oo Llinite) slid block m ‘ 
iF roe ( 
a= D CYA) aye), 
1s lYen b : 
J . ( Fyne) | 
| rw 


x4 = D Yen, My Yen) a oad al " = A 


wl Ly 2 an INnNnDdDow US QNtL Symécs 
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— WAPPialle 


pee T e 


Cope CoWSThCTION IDEA: Exploit sini lar stucure of FSTD and FSM 
ACS Bint 


E Kame ; Kote |:| Code Frown U5 she ona 


FSTD 
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O Beit tate 
ex eo, a 
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Envdey FSM | Decoder 
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be Gs ian Lad} Pave, ie a 
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More SPEQiFic TECWIQUE ; Rite YM: WwW Code how. U.S 
Find Graph BH owith fs owing proper Ties: 
) Each vertex has A out going edges 


2) Data LABELS? Each mm blocle - 


Appears AS a Iabe) a one edge 
From tach vertex 


3) Cope LABELS: Earh edly ¢ has an m-blotk Aow 
EE SS. The TD given by H with hese 
labels represents S 4 or a subsystem 


of (~ 
4) Deécopadyrity : 
) ates aa For k large enough , al| paths 


ONY ONE path SYMbee Ck — e., C, ¢; — Cp whith generte 
the same code Sequence have the 
Same data label on edge @, 
Eneenines Hck inrhal etate  Felled edges according te date 
mpats , reading of F code labels ii 
DECopnG: Sliding block , with lok- ahead and look-back of k symbols ans fur 


Defn: The L tate -transition matrix T corresponding te G is Given by 


= Onn ;* (Ff there are te tyes fom ete i 
t O i vo tdaes i salle? 
ADS ACeNcy WAT eX | 
The wwe sums of T Indicate th number of 20 coe 
owt Gein, eda es fomthe states of 6- 


Romar ki: 


Kemark 2: The state - transite matrix oF G" is T™ Riot. 


ee? 1o f Ye Mater WL Heals ; ae 
- oF oe aie Phohekly, The €06€ 


ENOUGH €06ES 


Sums sams 
> I< 2 code possible fiw) 
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Cope Consravcrnd Eeamme: Reto Id cole fom U4 to RLL(o 1) 
(F requenty Modwastionw » Fi a) 


\ ~~, \ 
C ~ aa 10 
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= Dis cond. extra age AL stale (> QV > > 
SSV0) ta lebels ‘fu | 7) i o , 


o/ lo . 
RK 
“ona on 
f i red ae i a o ‘i ep, 
SIO, CENCE RTEO WTA of 


\\ i 21; | ‘ TED WTA suet Cages 
~ Mae ‘nits © wh O te get stte-free enceder (block tode ) 


| ll 
O-1)0 
= Decoder i ob viovsly a block decoder 


Kemet : Code. Resign Involvec Chorces , discarding Cager 
ASMining data labele 
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C ope Rares AND Shawnon Capacity 


Dew. For a DNC the Shannon Capacity Ci deni 
bythe hmjit ! 
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whew |! 
# DATA SeOvences 
Aine of nen of | est wm in S he ws 
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i on, f 5? i din 
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Rework: Units are 5) its /oy kt 
Also use the os 09 Cap () 


Remark 2: Capacity Is also called. topological entropy (symbol it dy vam CS) 
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Copiwe THEOREMS 


Theorern (Shannen) LF there Gist a code at rate mm fom 
Un te S thee 


(vv) L< ¢ CAOACITY Is) UurimMate 
— 
Ww — Cap ( a oe Limit 


Theorem 2 f s hennon) For yates mm Von < Cap és ? Converse 
there exists © code fiom U, te S 
with o (finite) chiding block. becoder 
with rate mp: mh, Sor Some £ > | 


Code MAY NED TO G (au SYM Rot$ 
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CALCULATING Capac iy, 
Facrt: © re present the growth heke of the number of 
Stangmces of length m nv, 


Focr a: THE (+ 1s irre dang ble (ever date v reaches every state j,) 

Shannon 
1} One Shannen ‘Pamence. If generated. by more 
hon than one path Fem a given sate) 


C= Ay 
Where 


—_ GB 1¢4 
— ZF] 


Caview oF Lineage. A L¢eera 
ae —— SS ieee 


T-AT)v =e ae eae 
y ) ies We Bai Tey \ 
YET - ee Gar oer > any per\ G 33 

G3, G32 "73 4.2 G3 
-G,z, det Fe a) 


Requires A)et (T-AT) ile 
| ee ae etl +43, Det e pe 3) | — 


chor act tr. J stic | poly nom ; al PE bbol- Cbogenvs Thence 
a fife Ne 
AL (a1) T= [i ol vet 
—(-Aed)-| 
A=A- 
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Example 
dt (T- AL) = 
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\\ 
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fA 
CxPac iT o¢ Vaeinbe Lenomt GLAPHes 


Facr ae T + f S has Vari able length Synbolk, with 


[ Shannon \ ic et L. bis) Athen 


~ “Capl9) = = hog X, 
where X= Jargest real _* 
det has ) ia 
Enample : ELL (a, k) 


det 
00..01 Ata 
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: RIS AND BLOOMBERG: CHARGE-CONSTRAINED RUN-LENGTH LIMITED CODES 3455 


| TABLE I 
- CHANNEL CAPACITIES OF CCRLL Copes 
¢ 
4 § 6 7 gg 9 bb hl wp bb OO 


So RL 8 


5000 .6358 .6662 6778 6834 6866 6885 6898 .6907 6914 6919 6922 692 6942 

76648244 8468) 8578) = 8640 = 8678 = 87048722, 8734 874487518757 8792 
7925 8704 9012) 9165) 9252, 9306 = 9342, 93679386) .9399 94109418 9468 
7 $e 8832 9120 «9380 §=.9486 )=— 9552, 9596) .9627 = 9650 9667) = 69680 96909782 
ae * 8858 - 9256 9460 .9578 9652 9702 .9738 .9763  .9783 9798 9810 9881 


3471-63822 3931 -«.3978)»= «4003'S «64018 )~3=— 4027, 4034S 4038)=— 4041S 4044-4046) 4057 © 
4248 ($000 .§237 .$341 5396 5428 5449 5463 .5473 5480 5486 = 490.5518 
- -  §391 5746 © .8905 5989 = 6039S -«6702s.6093. Ss «6109-6121 = 6129S 6136) .6175 
= - 5497 5947 6153 6263 =.6328 ~— 6371 «6400S 642164366448 = 6457 6509 
.. -- - 6020 6260 6391 =.6470 = 6522 6557 = .6582 6602 3=— 6615 6626 §=— 6690 
-- - - 6039 6305 «6451 = 6540 = 65996639 6668 «= 66895708718. HZOD 
- -. -- -- 6321-6577. .6574 6638 6682 6713-6737 6753 —«6769 ~~ 6853 
. -- -- -- 6325 6488 «= 6590S «66057 6704-6738 «=.6763 = 6783 = 6798 «6888 
- - . - 6492689766666 6715 6751-6777 6798 =~ 6814 6909 
- -- -- - - (6493 6600 6671 .6721 6758 .6785 6806 .6823 .6922 
6601 .6673 .6724 6761 .6789 6811 6828  .6930 
3§3 2859 «2863-2866 = 2868 «= 2869-2878 

-- 347]. «3777-3893. «3950» 3981S 64001-64013.) 64022 .4029Ss 4034S 4038) 40874 
-- 3723-4199 4384. «447552604557 4578) = 64593 4603) 4611 = 46174650 
- - 4366 461447374807) 48514879 4899) 49144925) 49334979 
4418 4718 4870 4956 $011 .5047 5072 S091 = S105 S115 $174: 
-- 4761 .4935 5036 =.5099 5142 $172, 51945210 .§223. 3293 
-- 477449655077 81485196 823052855274 52885369 
. -- - - 4977, 097, $174. $227 5264 5291-8312 5328) $418 
- - - - s98) 5107 “SYBR™ «52445283 531353355352 .S.480 
-- - - - - - S110 -.5195) 5253. 8295. 5325) $352 53698471 
- : - - -- - S11 5198 = -.5258 = 5301) 533353575376 5485 
ot -- 1903 2101-2162) 2188 = .2202, 2220) 2228 © .2219) 222222232228) 2232 
| - 2434-2902, .3049 31123146) 3166) 31783187) .3193 31973200) .3.218 
yoo -- -- 3224 3464-3570 36253658 36793694 370437113716 3746 
aes. 33329 3660 «3807S 3885) = 63932 3962-3982. 3996) 4007S 4015 ~—«4057 
-- 374) 3929 4029-4088 )= 64127 4153. 41724185) 41964251 
3774 «63990» .4107, 4179s 4224 4257) 4279-4296 )=— 4309-4376 
- - 4017. 4149, 42304283 -- «4320 4346-4366 43814460 
- - 4025. 4170-4259 4318 4359 4388) 4410 44274516 
-- -- -- 4179 4274 «4338 «=. 4382 4414 44384456 4556 


Roe pe he rr hr ODDO OOOO °O 


4 
* 
: 


< os = >> 4182 4282) 4349. 4396) 4430-4455 44754883 
1278 1662, .1747) 17791794 1803) 1808 =. 181218141816 1817 1823 
ae 2271) 2480) 2559) 2597) 2618) 26312639 2644 20AT 2 6A9 2669 


AANA DAWNAUAUMAUAN SPE D HA BO Bw Ww WW WW le le NN PI tN 


aad8. 52522* 295): 3019 2305S: 3078" 3092 3103) - 2310 = F1Nd:. 3182 

a ee = st 29753162 3254 33906 «9938 3360) 23374 «=.3385 3993: AN 
to * ae a 3030) .3267) = 3386-3453 .3496 3523) 3543) 3557 3568 3620 
bes 2 = % ne 3316 0.3458 = 3840) .3592, 3626 = .3650 = 3667) 36813746 

| dT, ose ae is = 6s 33336 3496-3591 3650 = 3690) 3719-3739) 37553833 
a a an ns on : 3514 3619) .3686= 3731) 63763) 3786 3804 3894 
: 13) a ie ke 3520 = 3633) 3706) 3756) 37913817) 3837) 3937 
1 Gi ADT -sI4Si- 1693; AST... oT52h. 1827 51830-1595. SSA 30556, «2582 
7\ 17132054) 2160) 2206 2.2230) 22442252 2287. 2260) 2262) .2281 
i 8! = 2318 2499) 2578 = 2620 - 2644 = .2660) 2670) 2676 = 2680 = .2 709 
ae 2615 2672) 2786 28472883) 2906) .2922, 293329412979 
| 10! ce 2755 2903) 2983) 3030-3061) .3081 = 3096 = 31073158 
ii 2786 2965. 3063-23121 «3189 3185) 3204 32173282 
112] ae 996> 3109 BEB «3222: 3253... 3275-3292 - <3309 
13} ne ag 7 3007) 3134 3212, 3263) 3298) 3323 33423432 
7 || 09340 1217) 127913031314 132 13240 .1327) 1328 0.1329 1335 

8 |} : #8 1691 185019101939 .1954 19621967 = .1970 19721993 

9 |! : AOG63 - 2133 32237 2287 12318-2332: 3.2342 2349 2333 238 
10) : oe 2268 2418 2492) 25332589 2576 2587) 2594 2633 
1| : 2320) 2516 2614 2669) = .2703) 2726) 27412750 2804 


* Ordinary RLL values (no charge constraint) are in right column. 


ACKNOWLEDGMENT [2] A. M. Patel, “‘Zero-modulation encoding in magnetic recording,” 
IBM J. Res. Dev., vol. 19, pp. 366-378, July 1975. 

P 3] A. M. Patel, ‘‘Charge-constrained byte-oriented (0, 3) code,” JBM 

ihe authors wish to thank Pat Kocsis for help in preparation [ Tech. Disclosure Bull., vol. 19, pp. 2715-2724, Dee. 1976. 

the manuscript. [4] J. Hong and D. L. Ostapko, “Codes for self-clocking, ac-coupled 
transmission: aspects of synthesis and analysis,” JBM J. Res. Dev., 
vol. 19, pp. 358-365, July 1975. 

[5] T. M. Chien, “Upper bound on the efficiency of dc-constrained 
codes,” Bell Syst. Tech. J., vol. 49, pp. 2267-2287, Nov. 1970. 

C. E. Shannon, “A mathematical theory of communication,” Bell [6] P. A. Franaszek, ‘‘Sequence-state methods for run-length-limited 

Syst. Tech. J., vol. 27, pp. 379-423, July 1948. coding,” JBM J. Res. Dev., vol. 14, pp. 376-383, July 1970. 


REFERENCES 


he fe pana caput iets 
af s ati 


Comment” on units : 
| = Aog | py 
LF b= 2, Units ave binary abn, ts//tymbe 
or bits /symbe) 


Te bes, units ore Tmary oly its tym bo) 


etc... p 
© 


rs 


STanstiot Fropanies oF DNC 


Questions we might ask - 
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Table 9-4 lists several group and block codes. They represent a 
complete re-coding of the original code, and some of the rules for the 
re-coding are rather complex. The coding is in essence done by breaking 
the original code up in groups of n bits and mapping them into m bits. 

The pattern of the m bits can be predetermined and stored ina library, 
with a pattern for each combination of the incoming n bits. Such an ar- 
rangement is made for the popular 4/5 code or GCR code, and the conver- 
sion table is shown in Fig. 9-30. This code is one of many /BM originated 
codes that paved the way for today’s high packing densities of 6, 250 bits per 


inch, and higher. 


Basically, GCR uses the NRZI format for “1s and “0’s, but a restric- 
tion is added: there can be no more than two “0’’s in sequence (k=2, d=0). 
This quarantees that flux changes occur at least once in every three bit 
cells, and the variable-frequency clock need only be able to lock onto three 
pulses, corresponding to a succession of “1’’s, alternate “1”s and “O’s, and 
a “1” followed by two “0’s. 

The GCR-code is also advantageous when error detection and correc- 
tion is considered. The reader is referred to Ringkjb’s paper for details. 
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